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ABSTRACT 

PART L 

ine exacusnr en semnperature series expansions for 
the [sing mode Wispecitic heal are studied using series 
analysis techniques. The estimates for the critical 
exponent. 0 (Prem tae series on the Tace-center cubic, 
DOGV=Cenver= CubLC. 2nd. sample cubic Lacvtices are Tound 
DO MCOnVerce ands cme wanalysisy indicates thar, the critical 
exponent of the Ising model specific heat is a = 0.114. 
fais value ls seed isartreemneny with Utne present accepved 
value of a = 1/8. 

the exact: hich temperature Series expansion Ler 
the 2 model wspecii vem heavy oOnmiwwhe face—cenver cubic 
lattlee 1s analyzed. “The-anealysis olf “vuheespeciiac heat 
Series and) 2us derivative 2c ound vo be consistent with 
eeSpecirve Viecave index O10) — U0, (corresponding Co a loga— 
Titanic Sinculaci oy. sel er Une uloOnal horm Tor chee xy moder 
Sepeciiic heat. women is concisteny with the analysis, 16 
Drepelued. 

The estimaves of the critical exponent for the 
specitic heats of the Tsing and XY models are compared 
with experiment. A good agreement with experimental 


systems is found for both models. 
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PART IL 

A Mew vest ol oscaling theory in the critical 
région is proposed. The low. temperature series 
expansions foretne Two and three dimensional Ising 
model mapnetizarvion on critical paths or the form 


p 
Ce tee al 


iS studied Usilig Series analysis 
vTechnLoQues and the estimates of the critical exponent 
On” nese paths 15 comparcU with tie predilections, of 


scaling theory. "A @o0d apreement with sealing theory 


ie Tound in both two and three dimensions. 
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SPECIFIC HEAT ‘SERIES 


CHAPTER 1 


HIGH TEMPERATURE EXPANSION OF THE ISING MODEL 


SPECIFIC HEAT FUNCTION 


in this chapter a brief review of Some of the 
standard techniques used to derive high temperature 
series expansions is presented. A thorough review of 
the various methods of deriving high temperature series 
expansions for the Ising model has been given by Domb 
(E960). 0 Ts =nmot the intention of the author (6. cive 
a lengthy review of the concepts of graph theory, 
which is widely used in deriving the high temperature 
SGGIes expansions. (For a preaver insient. into such 
concepts, the reader is referred to Domb (1960) or 
Sykes, Essam, Heap, and Hiley (1966). 

ihe Hamslvenien Ol vhesspim one-—noatiel eine mode: 
WitMeneares t= Nel enborainiveracc 1 ons May em Widiys elimi 


the form (Domb 1960) 


N 
a) awe 0505 Sage O5 (aie) 
a Wes see i=l 
where ees Es tf) pseone span  vardableovassociated evith tie 
ee site on the lattice, m is the magnetic moment, H 


the external field, J 1s the interaction energy between 
neighboring sites, and N is the number of sites on the 
lattice. Thesotvariables take the values = according 
to whether the magnetic moment is parallel or anti- 
Daralled to the magnetic field: 

The thermodynamics of the Ising model are com- 


puted from the partition function, which from the 
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Hamiltonian (1.1) takes the form 


N 
oa a on) SPAS AU eter ete dala A oY for al 
q o,=t1,..-0,=#1 <i,j> 7? 1=1 * 


Cine) 
where the outermost sum is over the ia possible values 
or oO, for the N dattiee sites. 

since the o variables commute (1.2) can be writ— 
Ten 2s asproducs 
N 


Gas lgme SD. Cecoms @m ils CxO Uluteey) (so 
Si aks 2 Schl : 


where K =) J/kT and L-= mH/kT.. The 64.05 Satviei ve pne 


relations 
2 4 3 5 
Ves = = = e A = a A = good 
(o505) (o505) Sees (o50,) (o,05) (o50,) 
(aye) 
and hence 
exp(Ko,0,) = cosh K + O59; splgion 1 Cie) 


The flcsuepreauce site ln) (calle bes cxpandedsas 


follows (van der Waerden 1941) 


1 ©Ceoshn K«+ O40, sinh i) ewes = K ft (i + 59, tanh K) 


<1,J> <ip 
(a 
= flees peyeel [1 + (tanh K) ) 495 
<i,j> 
é 1.6 
+ (tanh K) ) O595,9,0)4.-J (1.6) 
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where p(r) denotes the total number of ways of embedding 
in the lattice all graphs of r edges whose vertices are 
alt of even degree. The pir) are in general polynomials 
Or decrec: in N. Any such graph gives a contribution 
of 2 for each vertex when the sum over the appropriate 
oO; is done. Hence the aN in (1.8) arises when one per- 
POrms) Che “Sum over al lestates.. 

Instead of (1.8) what one is really interested in 
ie vane dimensionless, Helmnoltz Tree energy per site 


defined as 
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+ lim{logf[l + pir vy J7N (2270) 
N>o r=1 
where v = tanh K. It has been shown by Domb (GOGO) ieee eu 
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Tne zero siield specitic heat per site voan be 
obtained from (1.11) by differentiating twice with 


FESDECU LO The “Vemperature, ine, 


2 
9 
C,,/kT = ae (-F/k) : (Calgakes) 


CHAPTER 2 


SERIES ANALYSIS TECHNIQUES 


2,.L “Critical rxponents 


In theoretical and experimental results one usually 
assumes that the thermodynamic variables have a simple 
power law behavior near the critical point of a ferro— 
magnetic system (or analogous systems). Therefore near 
the critical point it is assumed that the thermodynamic 
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where Xy its the ceitical point, y the critical exponent 
and A 1s) the critical ‘amp iicude, 

Ap DPeCa seman, Ceneral Cerinitlon sol, a cra vaca l 
expoMmeny toe describe the behavior near tne criuledi 
POImuueor @ general function f(x) is given by 
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where f(x) and Se) are positive. Of course the exis- 
tence of the exponent y does not mean f(x) is simply 
proportional to CBE One must always expect correc- 
pions Lermeanoren loner order. simce: (2.1) represents only 
Taenmdominany asympcoreie mehavior. 

Por a compieve. definitions ofall criti caitexponents, 
ineluding’ thoser used) by the; author; the readergisercter-— 


red to Fisher (1967). 
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Pee favo Method 


The ratio method for analyzing series expansions 
basspeen used irs i by Donb and Sykes (195742. end b) td 
6stimeate critical points and exponents. The method 
has been reviewed by Fisher (1967) and more recently 
by Gaunt and Guttmann (1973). 

tf the dominens singusarity of al function occurs 
at some value X endeis ol Gie form 2, erie tie 
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tend to 


uy, = te ~ a - Ey, (253) 


where ee Ek ae it the successive values, of u, are 
blotted against 1/n, the Location of the singularity 
Can be estimated by extrapolating bo the intercept 

m= =. The value of y Can be estimated from the Jimil— 


ine slope sor the plog.) If Xo has been determined by 


this or any other means a sequence y(n) of estimates 
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To aid the extrapolation of the sequence y(n) one often 
plots these values against 1/n and extrapolates to the 
intercept n = © to get an estimate of 7. 

in eeneral the series expansion for fl.) exni bats 
nou-—physical singularities: as well as: piysical ones. 
The physical singularity 2s usually the closest sinzu— 
larity to the origin on the positive real axis. When 
Pnere ore No non-piysical singular@mii1es mearer 1 oO sue 
Origen trhan them onvsical singularicy  cvhie ratios sare all 
POsleave and when the nen vhiysical singularitics tie 
Weil nourtside whe @erre le or “converpence derinedsby ve 
paysotcal singularituy, therm the ravios will converge 
Papltadiy, If Uneradoninant singularivy Lies on che negza= 
Goiverreal axisetne sims, Of ihe Series coehiicrents wilt 
alternate and when the dominant singularities are complex, 
tne slens= will bewirregular. For "these cases, ditterens 
techniques must be employed. It should also be mentioned 


that more complex singularities may arise such as 
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The ratio method has proved most useful when 
applied to high temperature Ising model susceptibility 
series, since in<these cases the dominant. singularity, 
1SPOn the positive real axis. The low temperature 
series presented in Chapter 4 does not in seneral fall 
into this category, due to the presence of complex 
bairs of singularities symmetric about the negative 
nea! ats and, Nearer The origin than the pivseicals cingn— 
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Jos peLYranstormacioNs Of ExXpens on Varilepics 


Late Ossi ble | iowever, in some: cCaosee CO or anci orm 
The physical singularity nearer to the oriein than) tne 
complex singularities. This is accomplished by trans— 
forming the series expansion variable. A conformal 


transformation of the form 


t 
E> > Geron) 


is®the type of transtormation most orven used. 

Wortis (1969) applied a transtormation of this form 
GOPtne anti ferromagnetic suscepuipility: series of, the 
ising. moael On the isec.c. lateice.§ Afver thus traus— 
forming the Curie point to infinity the resulting series 
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More general non conformal transformations can 
also De used." Betts, Elliovyl*and Diezian (107i) 
introduced a non conformal transformation of the form 

ax! 
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and successfully analyzed the triangular lattice 
Pec uay TONS eserioo with 1G.) In Chapoer waste 
Pesuits of The use (Ol several transformations of ‘this 
LoEemrarevdis cussed. 

A-ORanCheDeInL. Sanaularicy ot sube storm aCe...) 
Gan be jtranstormed imvo <avsample pole by aking the 
logarithmic derivative of the series for f(x). 

the form becomes 
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and aeratiovanalysis on (2,0) sean be used as .a. test 

Of show welll the series Ls represented by the form 

eo thieestnee frome(243) shewratviossshovid approach 

the constant value ee ii een ees lame te som x is avail- 
able a further extension of the transformation (2.8) 


of the form 


Goren eec) + log f(x) (2.9) 


teipossibles Ifer(x)ehas thes form (2.1) then F(x) 
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Should converge to y. If the series expansion of F(x) 
is sufficiently regular an upper or lower confidence 
limit can be put on the estimate of the exponent. 

Tavs can be done by evaluating the truncated series 

rom Cx at the erivical DOLEG Xoo using each successive 
CCetilervent oF the expansion to P(x). If the series 
ie regular, the mies value or B(x) can be used as an 
upper or lower confidence limit, depending on whether 
the sequence of estimates to eo) appear to be con- 
verging on y from above or below. This new technique 
iieserles aati ts) DEIN. Diy ehORWard Dy cnet der 
ance iS Sse an rae sanalysis Ob Severa = Series aia 


Chapter 4. 


2.4 Padé Approximant Method 


The Padé approximant method, which was first 
applied to the Ising model by Baker (1961), approxi- 
Mmeves a cUNnCtELOn by one Tarlo of oWwo polynomiake. 
Following the convention of Fisher (1967), the [2,Mi 


Page a0preximant vo a power series sis defined by 
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C[L,M] = shbetie nosth ie or thbiie (2.0) 
On ‘3 Le Ope ata Anco Caen’ PP 
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(ne Coeciricients Sogn Sepia Pos Ayo +++ Ay are calcu- 


lated by requiring the expansion of [L,M] to agree 


exactly with the given power series, f(x),up to the 
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order (itM) = Rh. waere Ris the order of the term at 
Which OGx) aie truncated. 

Tne Convergence "or Padé approximants. 1o. (> ) 
has been studied by Baker, Gammel, and Wills (1961) 
and Baker (1965, 1960). This problem will not be 
Cons dered bere but (uheir result, that yihes in. Nal 
diagonal approximants may be expected to converge to 
f(x), will be used. Essam and Pisher (1963) also find 
that the [N,N+1l] and [N,N-1] approximants seem to 
converge just as well and have used them to estimate 
Ghe- function. “The mature of Padé approximant is such 
that ue Can represen, a simp 6s pole wi a. 1 Uneuon 
SkLaculy. SOltus cca rallew tT poss rem tGnv laite wim 
thee 1uUNCcti1.On being examined so Une Singularity as or 
Chace norm. “Convercence Va one reson Ol esuchma, polewis 
rapid whilepin the regionw of a2 bDrancm point 26 as 
eons i cderab iy ys lower. 

In order to use the Padé approximant method most 
eGtectively for funceions Of (the form (2. 1)5 tne serves 
should be transformed Using (2.9), 4 process which 
eonverts the singularity into a simple pole. 

bE san accurave estimate of the crivical exponent 
Ve tsuavarlable., estimates OL the critical poine can be 


obtained by forming Padé approximants to 


a Gy) ie as a) es (2,11) 
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It -an eeccurate estimate of the writical point is 
available, the Series can be ‘transformed using (2.9) 
and the Padé approximants to F(x) can be evaluated at 
Cue Verivi cal pcan x CO Obtain estimates of the Criv.-— 
Cal index y. 

One can also use the Pade approximany method Uc 
calculate estimates of the amplitude. Taking Padé 
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CHAPTER 3 


ANALYSIS OF ESING MODEL SPECIFIC HEAT SERIES 
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As noted in Chapter 1 series expansions of the 
Specilicehnéeatein zero field can be obtained as a 
series in ascending powers ofeK’=J/kTorty = tanner, 
Me sextrapolationstof the criticale indices were 
found to be quite insensitive to a choice between 
these variables and therefore analysis is presented 
in the natural variable v only. The critical behavior 
obuthesisine model specific heat expansions has been 
studped on thesface-centered cubice sbody—centered cubic. 
simple cubic (Sykes et al 1972b), hydrogen peroxide, 
hypertriangular (Leu, Betts and Elliott 1969), diamond 
(Essam and Sykes 1963), crystobalite (Gibberd 1970), 
and tne octrahedral (Oitmaa and ElIioGt_ 1970) Lattices. 
SLACe UMC «cc. LavuLcen tos tne, MOosvecloce Ly 
packed of the lattices considered, more "information" 
SbCuUt Une, SUrUGCUUre 15 “Convaincd, in Vie vearlier ico= 
S2picients of the t.c.c, expansiom, | Hence 1) ss expec— 
wece Gaia bliecs sericecet orm tne f.¢.c.. lavelee wilieoapproacn 
Loco lami rine Delay tuormanOnme: TrapLaly tial Une mselr Cos 1 .Or 
the other lattices. In the analysis which follows, this 
abhi ceawil be tne (oe principally considered jasince 
the emphasis is on methods of analysis. However, most of 
the calculavions Mave been repeated for the ouner Javeices, 
and where the results are sufficiently good for infe- 


rences to be made about the behavior of the series. They 
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Support all conclusions based on the analysis of the 
TTC. Cs Lapuloe. 

The ratio method has been discussed in detail 
In OCCL On ge. ae bie CASS 01, ne speci icaness 


the assumed limiting behavior is of the form 
x — 16) 
C/K A Wy) : eo) 


Equation (2.3) for the limitine pehavior of the ratios 


Of successive coefficients of the series becomes 


a-1 -l 
Uy, ~ (1 + aa) Me 5 (372) 


Lf, ie LSU Eeedeahlwoctilavics Of (Cale bemouratned 1 om 
the slope of the ratio curve (when plotted against l1/n), 
and from the sequence y(n) defined by equation (2.4). 

To form this sequence the most receny estimate for Wee 
based om the hich temperavure suscepvulpilivy expansion 


(Synecner al Loyj2a)yussused. This esvimave is 


Vo = OTOL A0r0, 000005 


or Gs aap 


— 9.8290 + 0.0005 


The sequence, together with the sequence of ratios 
u, is given in Table 3.1. The sequence y(n) in Table 
2.1 seems to be converging to a constant value; for 


res Oo, ,)14 ‘the=valnes of y(n) are very close 
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R@vligg Or une f,C. 6. Spec ulac Neat, vorteuner 


with sequence of estimates y(n) of the critical ex- 


ponent usine the critical point il O LOIea 


3 8.0000 Ora 

4 8.2083 -0.3405 

Seuss 0.1935 

6 8.3458 -0.0946 

wh 8.5495 -0.0888 

8 iil -0.1113 

9 gee 7am ~0.1235 

10 8.9709 SoA 20 
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twogetner. The sequence y(n) Uselid dine) is pletced 

Ne cee 7 a Migure 250.) hlsowshown are eimiiar 
sequences for the bic .c. (dotted line) and s.c¢. (oreken 
line) lattices, which were the only other lattices with 
Feguler ratios. | or locse packed lattices ; such as 

tine D.c.c. and sict lattices, only even terms are 
Present the thie spectre heavy expansion ero (2.4) muse 
Dewslignuly mMedvitved. — These moditied. form is 
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(Sykes et al 1972a) have been used 
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Simple cubic V We Souk e 0 0002 


: (325)) 
COs 5st 0 0010 


Body centered cubic ve 


There 16 a close Similarizy oan the behavior of 
the three sequences; they seem consistent with the 
View boei o 1Or tne Llsine Model is Gevermined only 
byetche dimensionality Or the Latrice andinoy Dy aus 
detailed structure. They are also consistent with 
the view that a is very close to 1/6. To obtain 4 
Wannerot (4) forwthe: face cenvercc cubic, of J, lab 
requires: 1/v.- roe COs Thucr tt une limit 1s exc. 


1/8 the sequences must pass through a minimum. 
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Pesplovved against 1/n in the usual way. In addition 
in this Figure we have shown two lines (dotted) which 
Davens Opes Corresponding, to @ = 1/3 anda a = 0 ond 
which have the intercept we = 9.3290. The dast ix 
points ef the ratio plot appear to Iie on the line 
cCourecponding vo a= 1/76.s. This is a eraphical expres-— 
SPOnnOM Unew tact that tne: last six Values ot 7 Um) are 
Memy IC lOsc Onl) Oculus rates in veuinitcine 
Pasion how unlikely iG would be for the ratios, atver 
bavane Seveled down vo such regular pDehayior Which polmos 
TO.ene Known elnvercepr, UO Andicave Subsequcnicly anovner 
SIMpLe Valle wh on OO. sinew ratio MevnodmsuroOne iy suet esis 
@eicevery close to am 1/8 on the fic,c. Lattice. — The 
evidence of tlieeb.c.c. @ndvs.c. also tends VOmsuppore 
Dass assump ulon. 
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eulvswere very anconelusave. in Table 372 the Locavion 
Of one physical pole and the residue for a teow Padé 
approximants to the logarithmic derivative of the 
specific heab on the f.c.c. lavtice are piven, |NoGice 
that the Padé approximants have not converged to the 


Critical polnt @s weil as tie ratio method and also 
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Table. 3.2 


Estimates of ie and a from the series expansion 
for Cu from Uhe location of the poles and the residue 


of Padé approximants to d/dv log (C,/v°). 


ee Singularity Residues 
Elbe hla 


ee Oelolte 0.4225 
[5,6] 0.10256 0.4290 
fone 0.10264 0.4359 
er ats GeO 264 0.4329 
[4,6] 0.10284 0.4520 
ioe 0.10215 0.3932 
[6,4] 0.10289 0.4575 
esr 0.10340 0.4978 
C2561 0.10356 0.5052 
C455] 0.10440 0.5526 
[5,4] 0.10744 0.5486 
[6,3] 0.10465 0.5651 
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the critical exponent is not even close to the estimate 
from the ratio method. Since there is a good estimate 
Or ae Pade approximants to the function (2.9) evalua 
ted at v= 0.10174 will give an estimate of a. These 
ave siisved it Table 3.5. Notice Chav bney give a = 
On5> an Stronzecontrast to the ratio result. There vs 
also a large spread in the values for both Tables. 


Hunter (1968) found that the mimic function 
-a 
C./R At = OO eee al av/v,} (356) 
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Chis (eOnjecturesis right 1G would texnleinowiy direcs 
Padé analysis of the specific heat gives such poor 
results. |) [Meer ornee., Dae mimic sunmcoon sue rest 6 1 ine 
second derivative of the specifie heat with respect. to 
v is the function to analyze with Padé approximant 
methods. The first derivative with respect to v should 
also give better results when analyzed. To test this 
conjecture Padé approximant analysis similar to that 
done on Cy is performed on (d/dv)C,, and (d/av) °C... 

Tie results are givenein Tables 3.4 and 3-5. In Table 
So iwehe Location of the physical simeularity and ¢cor— 
responding estimates of a for the Padé approximants to 
the logarithmic derivative of both (d/dv )C,, and 


(d/av) °C, are given. These Padé approximants give 
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Estimates of a from evaluating Padé approximants 


to (v-v,) (d/dv) log(C,/v"). 


Approximant 


Value 

CLM] 
ee (ea 0.3499 
Rael 0.3463 
We 0.3601 
(7,4) 0.3440 
[4,6] 0.3542 
[oad Chee 210) 
[6,4] 0.3575 
Ei aeid Or 35 2) 
[3500 0.5563 
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PeSuleSs CuUIve COMSISLENL With the ratio method. In 
Table 3.5 evaluations of the Padé approximants to 
Pure Clon 2.9) ors he (d/dv)C,, and CUcCa mer are 
piven. Aeainm Doth eive results consistent wich the 
ratio method. From these Tables it is seen that Padé 
amalyeis Or the first derivative of Ghe specific meat 
Seems FO Sive TNE Most Consistent results. “Thails was 
also are Por sphesovhner seven Tsing model Se specihauc 
Mee besericowonaly Zed. la0le s. 5 ci1Veswenscspimapen at 
Gd2= O.219 from the analysis) of the first derivative. 
Lilie temyvery Clone sbOs UNC erat LO resui. 

(Mie Woe yanela View that une speci ie neat on 
a three dimensional isine model of a ferromagnes 
diverces ac the e¢ritical vemperavure, irom above, 
inversely as am one eighth power is consistent with 
the analysis’ presented here, but the author believes 
the analysic somowe OG vo De slightly tess than Wo. 
However the evidence that a is close to 1/8 is strong. 

In Figure 3.1 the sequences for y(n) appear to 
De inesr in 17m for the tasv few terns. (ineardy ex— 
trapolating this sequence to n = » yields y(~)=0.114. 
Mipow leet) ooreemeny winn Tanle S052 Uhusta walues on 
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CHAPTER 4 


ANALYSIS OF THE XY MODEL HIGH TEMPERATURE 


SPECIPIC HEAT SERIES 
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This chapter is concerned with the analysis of 
exact high temperature series expansions of the Spin 
1/2 XY model of ferromagnetism or of a quantum darcice 
fluid. This model was originally introduced by 
Matsubara and Matsuda (1956) as a model of a quantum 
latudleesiluid.- he spin 1/2 XY model is of ereat 
theeretical interest as probably the simplest quantum 
mechanical many-body system (excluding "diagonal" 
models like the Ising model) and it is also of experi- 
mental interest as a model of an insulating ferromagnet 
Oreancleerromagnen . 

MNeEChods Lor derivaraous, ol Bnerexpanstome was: 
WeOLeWemel Ven mclen Durie nL eres legs reagel ts sre ormed 
BORBbevis prt llobtmand Lee (197.0). and Bets, (197s). 
Bieven coefficients in the specific heat series. for 
Pict. crw Ce oOMGnDmCh Chu oun ces Navembecnm, Ce vecw py 
Betts and his co-workers at the University of Alberva 
(Betts, Elliott and Lee 1969. 19/0 and Betts and: Lee 
LOGS) ep Only the wie. c 5 Lattice will be studied jmere 
Samce the series cn the Dec.c. lattice gave very, poor 
Pesillts) and mo, conclusions could be made abouv ie 
Critical. behavior of the specific heat. on this datvice. 

From a racilo: analysis of the fluctuation terres, 
Bets, BLIiott and shee GLO 70) seStamatecdemme icici ca 1 
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On Uhe tf sence ml aLtice. THIiSseValueror thes ecritical tem— 


Perature has been accepted as the “best™ estimate of 
whey Crigical point Tor the specific, heavy seriec. as 
Tope casesor the Msine model the XY model will oe 
assumed initaally to be of the form (2.1). 

Figure 4.) (contaime a standard ravio plot stor 
ipiems De Cullen Ol ule fC. Ce MLavtice.  lhearata] 
Ci successive Coctitecients ie Te plOtLea “versus l/n. 
Also shown are two lines (dashed) which have slopes 
corresponding to ao = 1/4 and a = 0 and which have the 
intercept ee => 2e The ratios seknib1 Go very 
fLrOnsmOosci Lavon eame wiery lipule can pe esta mated 
Spouy cheir damivine behavior. 

When the Padé approximants were applied 
CEme cul yarvo she sSoeeidac heav the results were very 
ineonclusive. In Table 4.1 estimates of the location 
Orne er io Lead spol Ka and) the exponent of from tie 
poles and residues respectively of Pade approximants 
Pemtune) loparivamicuderavative ol the f <c.c. speci ive 
heat are given. This Table gives the estimates ie 


O2225 and a = 0.6. This is in strong conilict with 
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PIGURE 43a 
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Table se 


Bstimaves of K, and a from Padé approximants 


Ponta 7a) log(C,/K*). 


APDrOximant 


[L,M] pineuler. oy Residue 
So CR fas O2001 
ee O23 2 Ono02 
aor oe dees oe 0.603 
P68 On2 304 Oe 5us 
25d Cm 35 2 Olea) S) 
paid Ome 305 OnMOAS 
[4503 Ge2s60 On 622 
hore e 0.2364 0.633 
2 EE Oe 2 -3.98 
Seat ed ORs wee O02 O00 
ee endl Op2Sh2 O MOGe 
Berea Q,2274 O2202 
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Ko = 0.celUn trom Che tiluctuation series. “Also gq. us 
expected to be much smaller than 0.6. When the Padé 
approximants to function (2.9) are evaluated at 

K, = 0.2210 to get an estimate of a, a value of a = 38 
ts obtained. This is clearly not a very good estimate, 
and it indicates that (2.1) may not be a very good 
Gnoice for the form of the specific heat. 

In an effort to improve the ratios several trans-= 
formations of the form (2.6) and (2.7) have been tried. 
No Euler Cransiormation ampreves the ratios but 
several transformations of the form (2.7) do improve 


the ratios. The three "best" transformations tried 


SOS, 


GA or aa aes Hes) 
Le geek! ee) 


ey ce sae eT Gir sy) 
ly Dieta ie 


! 
k= ——— (4H) 
tO ie) 
In Figure 4.2 the ratios corresponding to these 
transformations are plotted versus 1/n, In order to 
plot all three on the graph, the ratios are divided 


t 
by the transformed critical temperature Ka for each 


FIGURE 4.2 


RATIOS ue VS. 1/7n FOR THE TRANSFORMED 


f.c.c. SPECIFIC HEAT USING TRANSFPORMA= 


MONG. (2). 
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transformation. The circles correspond to a = 1/1.4, 


b 


! 
1 and Ko 0.2844, the triangles to a = 1/1.5, 


besa and ie OFS0i4s ana che x's) toa = 1712 0.5 


b = 1 and Ko 0.3179. Also shown are three lines 
(dotted) corresponding tord = 1/450 = 1/85 and a= 0s 
and which have the known intercept of unity. Notice 
net in these ratios the oscillation is still present 
DULL considerably damped. “A precise estimate of @ is 
ROt Posstole. buvvall three sets of ratios appear to 
OsCittave abouy 2 value of a. < 1/3. 

Padé approximants to the logarithmic derivative 
of (d/dak )C,, and (d/aK) °C, do not give a useful estimate 
Het oie se 2Opr On Malus CO MOu LOCaue  uiae ny sacad 
singularity with any consistency. The Padé approximants 
to Ghe Tune tiem (2.9) for (d/dK)C,, and (d/aK) °C, give 
esvuimaves Ol O Which are Glearly in error. 

Tne specaric heat Series for tne XY model scene 
imtractable to the Standard Padé approximant techniques. 
When the ratio plots in Figures (4.1) and (4.2) are 
re-examined & pOSSible explanation for this failure of 
the Padé approximant methods is found. The ratios in 
both Figures seem to oscillate about a = 0. Ifa = 0 
the assumed form (2.1) for the specific heat is wrong. 
Tier orm 
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would be ‘more appropriate. Tf this assumption is true 
tue Cirst’ temperature derivative of the specific heat 
would have a simple pole and Padé approximants to 
(d/aK )C,, Should converge rapidly in the neighbourhood 
or the critical point. In’ Table Wl2 the onysical root 
of the Padé approximants to (d/dK)C,, are listed. The 
approximants seem to be converging rapidly towards the 
EBON C TOMS Geil. Jeleniane i 0.2210. The "best" estimate 
of Ka PPeom vale Leable sais Soe O.22033) “This estamatre- is 
remarkapl yc Osewrro Ene estimate from the fluctuag.on 
Sec lee.. =) Ver Converrence an this. Map lewis extraordinary. 
especially when compared with the convergence of the 
esulimabes Of Ko in Table 4.1 which assumes a function 
Ote newt orm: Vee meer Mule crit ee lero om mac Ned ila 
(d/daK)C,, must be very close approximation to a simple 
Dole inthe cribacal reston.. and the ay model Genecuine¢ 
heat will be closely approximated by (4.5), except 
DOSSIOLY Vieryoencer the critical point. 

Df thervmorm (155) ts. assumed, then Cstimates of 
Che crivieal amplitude A can be obtained py evaluating 


Padé approximants to 


(K - K,) (d/dk) Cy (4.6) 


Buevne crippled! Colney eS O,2210. Estimatesmare iven 
in Table 4.3. A best estimate A from this Table is 


A = 0.254 + 0.06. Hence 
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Table 42 


Estimates of K, from Padé approximants to 


(d/dK)C,,. 

Sete ie Physical Troov 
[4,6] 0.2199 
oye! 0.2199 
Meee 0.2208 
Eiger 0.2205 
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Table 4,3 


Estimates of A from evaluating Padé approximants 


to (K - K ,) (d/dK)C,, at Le Crean Gr 


Approximant 
[L,M] Value 


OF 

7 
[5,4] Oe 555 
fosse 0.2546 
ona Oe 2542 
C4 ,4) Oe Se 
Een 0.2596 
6,20 0.2554 
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C/K ~ -0.254 log(1 - K/K, ) =O 9254 CS : Cher) 


A possible simple "mimic" function for the XY model 


Specie pear us 
Cy/k ~ -0.254log(1 - K/K,) aL K/K J + o(K) (4.8) 


where $(K) is the correction polynomial. 

The proposed functional form (4.7) for the limit- 
ne behavior of the XY model specifie heat is tvenvavive. 
More terms for the f.c.c. specific heat series will be 


needed to verify this assumed behavior. 
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#3 


De. Che-ising Nodel 


REGeEnvuLy (Huiskamp 19 Te) ood agreement has been 
found between experimental measurements of the SveCiia¢g 
heat on the Ising-like material EP. and numerical 
Predictions tor the simple Cubic Ising model specific 


HNeau. Ihe experimental results were found to fit 


rather well to the curve 


= 


Cee a7) hayes Scale 


WES Te 2) 1.14 K is the antiférromagnetic Neél point. 
ines is Che same funetional form as sucesested By 

Hunter (1968) for the three dimensional specific heat 
CoCl also 


3 D 
appeered to fio the Punetional form (5.1) very well 


Cer wines Isalaver jilocloinn Wels sieietesleiye. igershes eur CS 


(Huiskamp 1972). 

These are the only Ising-like materials which at 
present seem to give an estimate of a = 1/8. Most of 
the other Ising=like materials seem to be fitved better 
COea Curve Wildes 65) (Cooke eual 19/2 )_ here 
seems to be some experimental evidence for a@ = 17s bus 
at present the evidence is not overwhelmingly in favor 
Gia tiie value. Hopciubly Hiener resolution and bevver 
estimates of the critical temperatures in experiments 


Will soon resolve this problem. 
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4 
5.2 The XY Model 


Bee near the superfluid transition temperature 
Ty is expected to behave like an XY model system 
(Matsubara and Matsuda 1956). This substance has also 
Deen examined in more detail experimentally than any 
ofner “XY jJike™ material. In fluid systems like He’, 
the heat capacity at constant pressure is the analogue 
ofetne  speciiic heat at Constant fMeid in macnetiic 
systems. 

Ahlers (1969) has found that the experimental 
curve Lor Cy Goes nov disagree with ao = 0 and) hist vnrece 
QLeernavrLvVe AnLerpretvavaons O17 the dauasal 1s had sas 

4 


OeUCIy we Ore recent sexperinenval Lresulus (Onesie Dy: 


Ahlers (1972) (private communication) fit very well 


(assuming a = a') to 
_ A Hi -a 
Cos = fos br et = ae 
p a 
for ee 
O203 (08 = 00 = 1002 
and 


One soe < oO 


This is a very good agreement with the assumed 
theoretical behavior of the XY model. The results of 


Chapter 4 are not precise enough to ruULe Out, Une 
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possibility that a is a small power and experiments 

do not rule out that a = 0. Experimental and theore- 
tical calculations seem to generally agree, but. work 
meeds to be done in both fields to resolve whether o 


Lovossiell Dower Or 2S Equal to Zero. 
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A NEW TECHNIQUE IN THE ANALYSIS OF EXACT 


SERIES EXPANSIONS 2N LATTICE STATISTICS 


CHAPTER 6 


LOW TEMPERATURES SERIES EXPANSIONS FOR 


THE ISING MODEL OF A FERROMAGNET 
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In this chapter a new technique in the analysis 
of exact series expansion in lattice statistics will be 
presented. This technique is applied to the exact low 
temperature-high field expansion of the magnetization of. 
tae sing, model. 

The low temperature series expansion is an expan- 
Sion about the ordered state. For temperatures slightly 
above absolute Zero the ordered state will be perturbed 
by thermal excitations. The probability of any pertur- 
bation from. the- ordered state iS given by the appropriate 
POLE ZMann actor.) sine ceneral, overturning of almost any 
Spal Causes ane increase anvenercy and the mosy inportanc 
perturbation--at-the lowest temperature will -correspond 
Gouaerelatively “bew voverburned spins. Then one canpeceeup 
the perturbabrons conveniently according to the number 
OG Overturned Spins, the enersy of any particularsper— 
GupbatiLen depending on the relative positions of these 
Some More precicely, the increase in enersyeot a 


DerbUreabion is civen by 
Ne =O2U (gs = 26) emis C6s 1) 


where gq is the coordination number, H the applied field, 
s the number of overturned spins and r the number of 


nearest neighbour pairs in the overturned configuration. 
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Denoting exp (-23/k,T) byez) and exp (-2mH/k,T) by 
Hu, the Boltzmann factor corresponding to (6.1) will be 
ak aia ie At sufficiently low temperatures both z and 
uw will be small and the partition function and free 
energy can be expanded as a double series in powers of 
Zand Us 

Denoting Noas the number of ‘sites: on a letiice, 
the number of perturbations for the ordered state cor— 
responding to a given Boltzmann factor will be a polynomial 
Piel. Une cont lecuravional free <enercy per yspin ss epRepor— 
tonel vo che lOcarithm of Tie .“CcOnliecurationalipartatvion.,.unc— 
tion and it can be shown (Domb 1960) that this corresponds 
VOuLFaki ne Che =socthlelent OF tne Pirsy power cf ei Gamecne 
Povure TOM Cuneo LOnew Hence. JeETLineg vole, linear parusct, 
the total number of Ways! Of choosing S Spins with 2 “bonds 


Dymlcer!, then the logarithm ofr the partvivion=funcuten 


per site is given by 


GQs—cra.s (6.2) 


log A = ) Resa u 


ais. 


It 16 customary to eroup the expansion CGO eee 


series either in powers of np, written as 
log A= J Lg(2) uw (6.3) 
Ss 


where L_(z) are called the low temperature polynomials 
Ss 
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and are finite polynomials in 2, or in powers of Z, 


written as 
log A=) (nu) 2% (6.4) 
n 


where the vO) areniinive pol ynomialeringise Yemmeetiic 


nOvaG@onnuheyimecvenergey per spin Will then be siven by 


Fe-5qJ-mH-kTloga . Coney 


When using the method outlined above, one usually 
RrPoupse Une DervUuroavulons accordine Gor bne Lopolocy oF 
Ulcinm nearesy MelpuDor = binkagces. “Tor tililustmave. one 
Ceol eci some niOd.  CONSIader, Ge, Per uUrDav dono Ol soul 
spins and three nearest neighbor bonds on the f,c.c. 
PacGeice. —(4e3%seonsicus of The Tollowing pervurbatious 
together with the number of ways of putting each pertur- 
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Hence. the oerficient of zt? he te 126. Many 
very sophisticated techniques have been developed for 
counting such low-temperature confisurations. (Domb: 1960, 
Sykea et 21.2965. 1973). 

The tests in this chapter were made possible by 
recently extended data on the low temperature expansion 
of the free energy on the two and three dimensional 
Ising models (Sykes et al 1965, 1973). The low tempera- 
ture series expansion on the Gwo dimensional 14tcices 
are now complete to the L,, (2) polynomial in) tie Go.) 
grouping and the vig) polynomial in the (6.4)° grouping 
for the honeycomb, Ly 5(2) and Voo(v) for the square, and 
Lj (2) and ¥35(u) for the triangular. The three dimen- 
sional diamond lattice is complete to Lj 762) and Va9(h). 
Mhe tests were also used on the hydrogen peroxide lattice 
data, developed by Dr. D.D. Betts and his group at the 


University of Alberta. The hydrogen peroxide hac wae 


is complete to L5, (2) and Wig h)- 


| | ie 
pce l shins) saga ature 


aor aaah? 


$i 
es 


e 


eer ,eabs 
7 7 : Por 

ya Sidiwead “es POR eaten sti at. otaed peal 

’ = : 

golsurase: 5 sisamae Wak ety RO ey AD buties ee Clam 

fenatounorin #étal. te oat ashy co yah sie, oss | 

-Ptaqme? w: a CewepAmees ia: te sedge) rae 


noetitel )Lahereisett- omg On 10. fuera ee asl. 


ba 
a 
i] 
~~ 
.— 
>) 
7 
Pd 
al 
pa 
¥ 
— 
< 
< 
. 
i 
—_ 
= 
ba) 
® 


; . 7 
eORish jae ad’ ert reigns te) alt yl ha tel fae ts “ust 
Lamy 
baie) Chi adi oT Seles gs urligee ‘aete.hp ferro) 7 


eatsjal sURzoa0d (pee rUr, "aa Ae Ses ora otaet ont : 


_ 

Pa ql 
<) hj 
euiOje!l. sbirte ten meee ving ant e7720fk So tetenevtnll 
-t 94 adelqmge: # 


Ose QueTR) Bhd hae Ayead: 1676 oo bogotewen a 


CHAPTER 7 


SCALING AND THERMODYNAMIC RELATIONS TOR 


INDICES 


ens 


py 
{.1 Thermodynamic Inequalities 


The only rigorous relations thus far proposed among 
the critical-point exponents are a set of VHegitva Vi tases. 


Using the stability-convexity relations 


( ) °F 
C,(H,T) = - T (S35) 20 Cia) 
M yt? 
and 
2 
at jE 
= (>>) >0O Gia2) 
Xp (HT) yMo T 


with standard thermodynamic manipulations Rushbrooke (1963) 


rigorously proved for any system 
Caer rie. 4 ay ierse se. |. (ie ® 


(We use the widely accepted notation for critical expo- 
nents; see Fisher (1967), for example). 


A further important inequalivy 
se er (I Ges) ee 2 Giae) 


was rigorously established by Griffiths (1965), by the 
use of very general convexity properties of the free 


energy. 


For the Ising model, Buckingham and Gunton (1968) 


proved the inequalities 


pe 9 = d6e= 2) 708 2) rae) 


, ig ‘ Ca hie «1 bu 


‘ 
oy ig 
Lava with fo" caved gm 

“ = 

~ 

; wey er 

= sho Tr) 

Dt a 


wis we Peae nh) avitriignd wes haz Lf ee 


ay yatxeuric 5 fet * 


53 
and 


Bee Oye + ey Ce yiiG es tart) (7.6) 


where d is the dimensionality and more recently Mishier 


(1969) showed 
¥ s (2 = y)v : Clal) 


Very recently Gritfitnhs (19/2) proved form tic 
[sing model that the magnetization on any path, which 


Mies UnemCh vical rer On is tom che Porm, 


eo ee BE ; Gio” 


Must Nave a cravical. exponent Witin a value less thane or 
equal vor. INE Drool Of this anequadity LoLlows -direer— 
iV onmr lew COnvemt by Olemhe Gibbs pireecsenereay. 

ue M, (HLT) denotes the magnetization on the path 
defined spy (7:8) B, the critical exponent of M(H,T) 
near the critical point MOL, T) the spontaneous magne- 


tization, ana © the reduced vemperavure CD ee then 


B 
eg Ga) Sta P | Gis 
and 


M(O, ; 0 PF : (7.10) 


Since the Gibbs function is convex 
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fhis implies the projection of the slope in the H direc- 
tion, on the M(H,T) surface is always positive. There-— 
Tore along any isounerm 2n the posative fH direction the 
magnetization will remain unchanged or increase in 
Vewvlen | (iis eeshown seometrically in Mieure (li ae 
sOl10d curved line represents the spontaneous macenetiza— 
Gio, stie broken curved line as any patm om tne form 
Gio) and T, is the critical point. 


Therefore 


M,(H,T) >M(O,,7) . (7.12) 
RauacLonse (29). C7 el0).. and (fei je simp ly 


: CG1ay 


Hence, as tT approaches zero 


8 2 Bo janes seniuelah ee Cree) 
In a private communication Griffiths has stated 


that he has shown from the Kelly-Sherman inequalities 


(Kelly and Sherman 1968, Griffiths 1967) that 
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By making some quite plausible but more special 
and less fundamental assumptions, a variety of further 
inequalities can be derived (Stanley 1971, Stephenson 
1971). The above inequalities are of particular Lae 
when used with results derived from series analysis and 


experimenirs.. 


ioe soca ling Theory 


it has’ been conjectured by Essam and Fisher (1963) 
GiateGl.5) can beereplaced by Une equality bute so rar 
NO) Cieorous, DrOOl "exisus alvnougch the conjecuure is con— 
SLSvUehy Wilth Gxperimentbal and model calculations and une 
non=rizgorous scaling—-law vheory of exponents which predicts 
woaveuie Inedqualivies ((.so-(.ls (eles (al>) are equalivies. 

The scaling hypothesis was suggested by Widom 
OP965easb) Dombeena Huncer (1905), Kadanott et 2) (1967) 
and Patashinskii and Pokrovskii (1966). The basic postu- 
mace wor ne  stavic eocakane hypovnesis asserts (nau one 
Gibbs jpetential Giit.h) is a generalized homogeneous func= 
pieiew A Lune yon 1 .8..)) es Dy definition homogeneous af 


MOteat Lovaluces ot etne parameter A, 
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Thus from the general definition of (7.16), the static 


scaling hypothesis states 
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Diprerentia,inge both sides of (7.17) with respect to the 


Piele Cerivative HH, 


M(H )e= 12> MOCH. yes). (7.18) 


: ul 
Letting A = (1/7) Ae and evaluating (7.18) along the 


exis Hi= 0, 6 1s found to be 


32 Sse : (7.19) 


l/fa 


Seuulne =) H in (7.18) and evaluating the equation 


BvOnM at Ne scribe Cale-SOLMe UM lus 1S.) POMnCmr ia) & 


Sr Sas . (7.20) 


Bqauatioms (7419) and (/.21) scan be solved simultancousiy 
for the scaling parameters. a and b, and then substituted 


Mico w/.lo) to ClVe the Magneric equation sol statue, 


M(H tT) =e NOCH, 7 eo). Ce 
Vee 
eer neue) uew clidia ss 
, ed ® Bas 
Putting A} = T o@ } =i two alternate forms for 
Cine sare possibile, 
Gry 18 M, (H/t") Cieaee 
and 
M(H,t) = HY? My (t/H?/*) Gime 
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where My and M, ave analytic in tne vicinity of the 
origin and A = 8é. 

One can obtaan additional exponents by taking the 
second derivative of (7.17) with respect to the field 
to get the isothermal susceptibility and by differen— 
tlating twice with respect to temperature to obtain the 


Spetiibe heat as constant field, This yields 


Y= yo tee 1) (7.24) 

BY ke Se, Se Mee Cie) 

Ohareteee ty erro. ae CP.26) 
and 

fe=ao! ; pee) 


i@rshould be clear how to obtainvall Che ici vest 
Sxponegus im terms of the scaling paramevers a and ale, 
and how these two parameters are eliminated to obtain 
a whole host of equalities among the exponents. It is 
alsoeioted that the inequalities of section -&. sab 


become equalities in scaling theory. 


fo eelesus Of, -sea linge neory 
For experimental magnetic systems or, with appro- 
priate identification of variables, for fluids “the scal— 


ing hypothesis can be tested by fitting a2 scaled equation 
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of state to the data (M. Vicentini-Missoni et ine oOo 
and references therein). More sensitive tests can be 
achieved with theoretical models, in particular the 

Ising model. The critical exponents derived from series 
analysis techniques can be used to test the various 
equalities (equations 7.24-7.27 etc...) predicted by 
scaling theory (Kadanoff et al 1967). The scaled equa- 
tion of state can be constructed by the use of analytic 
continuation methods on series expansions and then the 
tneorecical expectations: resulting from the form of the 
equation of state are checked against series expansion 
results (Gaunt and Domb 1970). The critical behavior of 
the magnetization and its temperature derivatives can be 
examined Om Une Crivical tso0uherm and) the estimacves 

of the critical exponentS compared with scaling predic- 
tions (Betts and Filipow 1972). The series expansions 
ORebnes higher derivacvives of tne free energy wie respect 
to the magnetic field can be studied above and below the 
eritical point to verify scaling predictions (Essam and 
Hunter 1968). All these tests have given good agreement 
with scaling and none contradicts scaling. Relations 
among critical amplitudes are also obtainable from scal- 
ing theory (Watson 1969, Betts, Guttmann and Joyce 1971) 
and they too seem to be satisfied by modeis such as the 
Ising model (Gaunt and Domb 1970, Betts and Failipow Lose). 
The next chapter will put forward a new test of the 


predictions of scaling theory. 
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Equations (7.22) and (7.23) yield predictions of 
the critical exponents of the magnetization and its 
derivatives along the t and H axes and these are the 
exponents usually examined and tested. The above 
Scaling Gheory equations also tell us something about 
the magnetization along any path which in the critical 
FesLGnsico of thertorm (€7.3). 


Along This path Give2) and (7 «23 ) become 


MCHA) 


cM, (> PAI/P) Cena) 
and 


M(H,t) = Hl/® Pela ceaan (8.2) 


From equations (8.1) and (8.2) it can be seen thac 
PMemduacranvais>. Om elt 280 wis Uivided 2nvOn tWorrectole 


= ql /BS | Hormecurves wach. < l/ Boa 


by the curve T 
exponents have their t axis value as given by (7.22) 
While tor p > 1/766 all exponenvs have their H-axis value 
determined from (7.23). 

These predictions from scaling theory are much 
stronger than the thermodynamic inequalities (7.14) and 
(7215) by Griffiths. It is noted that thermodynamic 
inequalities have once again been shown to be equalities 
in scaling theory. The predictions can be tested in the 


case of the two and three dimensional Ising models where 


extensive data are available in the low temperature-high 
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field expansion of the free energy. From Chapter 6 one 
seesthat the free energy is expressed as a POWEr Series 
iieune variables: uw = exp (-2mH/k,T) and z= exp(-2J/k,T). 
ives convenient To replace = by = = 2/2. in the expan- 
sions. Figure 8.1 depicts schematically the (s,y) plane 
showing the critical point, (1.1), the path i = 1 ‘and 

Ss = 1 along which the behavior of thermodynamic functions 

is usually examined and the critical curve s=1- Gene are 
(broken line) dividing the area of interest into two regions 


and the paths 
= - # p 
S a (Gl u) Cons) 


Btoneawolche The Critical behavior of the magnetization 
Was investigated in this thesis, 

From (6,3 )eand ©6065) and the detinition of the 
magnetization, the magnetvizativon in the variables weand 


ZS ea Oy 


M(u,zZ) = m€1-2u(aL/au),J . GD) 


Using this definition a computer program was 
written to derive series for various paths defined by 
(8.3) from the low temperature polynomials of the honey— 
comb, square, triangular, hydrogen peroxide, and diamond 
lattices. In the cases of the square, triangular and 


2 
diamond lattices s was defined as s = u/u = (2/2) : 
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PATHS ON WHICH THE MAGNETIZATION WAS STUDIED 
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For, the derivation of these series Ges Ae aMnieda spe expanded 
by the’ binomial expansion. If p is non-integral, the 
truncated series expansion will be a Very poor represen— 
tation of (8.3) near the critical point. For this 
ncesonmonlvwentverral values of pecan De tised in u(G.c0, 
iis, widel Wee ve sseries: for the diaronal path s =u and 


above. Below the diagonal path (8.3) must be inverted, 


{SaaS cree yee (8.5) 


this jean be expanded for integral values of 1/p and: series 
beitow the diagonal path can be studied. A second consi— 
déeération is’ that for a curved path, the closer the patn 

is to the diagonal path the more of the expansion data 

ws pused. Forothis reason only small inteeral values, of 
prend i/p are studied. 

The paths along which the magnetization was studied 
for the five Lattices are shown schematically im figure 
Sabeand these «ane Ghe paths for p= 3); c2.. lyei/ey and 
iy/ 2 ep buswauhe-cripvcal isotherm and szero eile ld avis. 

For the two dimensional Ising model 86 = 1.875 and for 
the three dimensional Ising model £6 ~ 1.56. Thus, 
scaling theory predicts the magnetization on the se ehpaeay) 
isotherm and on the critical paths for p = 1, 2, and 3 


must have critical indices equal to 1/6 for both the two 


and three dimensional Ising models. The zero field 
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magnetization and the magnetization on the critical paths 
for p = 1/2 and 1/3 must have critical indices equal to 
petorsbDoth models. 

The critical points for the two dimensional lattices 
are known exactly, so all the series derived for these 
favtices are exact... Nor the three dithensional lattices 
Che critical points are not’ known exactly. The locations 
of the critical points have been determined by the 
aubhor from a reanalyeia of the high temperature suscep- 
Giles tivy series Of boun Lattices (Hssam end Sykes 2963. 
heuw, Setvus and Hiljett 2969). The critical pointe used 
were Ae = 0.5.0/393 for the hydrogen peroxide Jiaturece aad 
i at ee 0.227832 for the diamond lattice. 

The series derived by Chis method are very Long 
Onesone of thespabhs.. This Venctn 2s. deceptive since 
phere 16 actually less conficurational intormationm an 
these long series than there is in the series for the 
spontaneous magnetization or the magnetization on the 
eritical isotherm, which are much shorter. For instance, 
on the honeycomb lattice the spontaneous magnetization 


Leecompleve vO gi? the magnetization on the critical 


isotherm is complete to tone and the magnetization on 


the path for p = 2 is compleve vo 537, The magnetization 


series for all the paths on the five lattices are given 


in the Appendix. 
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FIGURE 9.1 


RATIOS Un VS. 1/n FOR THE MAGNETIZATION 
OF THE TRIANGULAR LATTICE ON THE CRITICAL 


ISOTHERM. 
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ratio method agreed with other methods on all series 
except the diagonal series. The reasons for this 
discrepancy will be discussed later. 

Four different techniques other than the ratio 
method were used to analyze the series to obtain 
estimates of the critical exponents and three of the 
tests gave very precise and consistent results and 
the fourth method gave better confidence limits, on 
BOMe Of whe series, than the ratio methods The three 
tests which gave very consistent results were all 
variations of Padé approximant techniques. 

The first method consists of the determining of 
the poles and residues of the Padé approximants (Baker 
1961) to the feogarivhnmace derivative Ol tne magnevizauson 
Oipeone sviarTous Pauls. “ne Poles solve cavlmaves=o1, cite 
erm trca We pOInb anu Une residucs ive Gsyllaves On suite 
eriserucal expenenc. “In Wable 9i2, the locatloniem the 
pole and the value of the residue for a few of the Padé 
Aporoximants are tabulated for the path s= 1- Gen 
on the square and hydrogen peroxide lattices. It was 
noted that the location of the pole and the value of 
the residue seem to follow a smooth relationship for 
all Padé approximants. When plotted the two estimates 
form a very smooth curve for both series. The inter- 
section of this curve with the line s = 1 is a "pest! 
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Padé approximants to (d/ds)log M(s) on the path 
L 
Ss = Z/z =l— (1 = )~ on the square and the hydrogen 


peroxide lattices. 
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Pecheoh=O18—-0.01 i= u)? on the triangular and hydrogen 

peroxide were very similar to the curve shown in 

Figure 9.4 which is the same plot for the honeycomb 

diagonal series. The intersection points are found 

to be 0.1265 on the square lattice and 0.319 on the 

Hydrogen peroxide lattice, The plot of location of 

the critical point versus the residue forms a very 

Peculiar eand smooth curve for ali the other serics and 

appears to give @ very precise estimate of the exponents. 
The second method consists of determining the Padé 


approximants to the series for 


CL) (adi). Loe eM (i) 
Oe 

(1 - s)(d/ds) log M(s) 
Oe 


(z - z2)(d/dz) log M(2) COLL) 


and evaluating the Padé approximants at the critical 
eis (GIS Ie SE Beko ay Zz.) tO give estimaces of 
the critical exponents. Table 9.3 represents typical 
Padé evaluation tables for the two and three dimensional 
lattices studied. From this Table a "best" estimate of 
the critical exponent for the triangular and diamond on 
the path s = 1- (i = ee ts) 0.12495 and G.31> respec 


tively. The scatter of the numerical values of the 
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evaluations of the Padé approximants to (9.1) in Table 
9.3 are typical of the scatter in the Padé evaluations 
ofall the series tested. It was also noted that the 
higher degree Padé evaluations appear vo: be converging 
towards the scaling predictions. In all cases the 
results from this method agree very well with method l. 

The third method involves making Padé approximants 
to various powers of the magnetization on a given path 
ang chen plotting the location of the pole versus the 
power for a few of the Padé approximants. The intersec- 
tion of these curves with the line corresponding to the 
Critical point gives a "best" estimate of the exponent. 
Pigure..9.2 is such a plov forrthe series on the patn 
s=l-(1- We for the square lattice. From this pilot 
a "best" estimate of the exponent is 0.0654. This method 
gives results very consistent with methods 1 and 2. 

These methods were the only methods which gave any 
Coneistent results for many of the series. 2 fourmsn 
method which worked successfully on the series for the 
zero field magnetization and the magnetization on both the 
critical isotherm and the diagonal paths was the numerical 
evaluation of the series to the function (9.1) at the 
Critical point using successively hicher Goch imGilenie 
of the series. This series should form a sequence of 
estimates converging on the critical Exponent. Ga vimt oe 


sequence is regular the last value can be taken as einer 
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Evaluation of Pade approximants to (l-s)(d/ds)log M 
a 
eu the Crities! pointe, s = 1 for the path s= Z/Z = vies 


on the triangular and diamond lattices. 
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an upper or lower bound. These sequences for the magne- 
tization on the critical isotherm for the honeycomb and 
hydrogen peroxide lattices are shown in Table 9.4. From 
these sequences it can be estimated that 1/6 > 0.0659 
for the honeycomb lattice and 1/6 > 0.1840 for the hydro- 
Ben peroxide Tattice. When the inaccuracies in “the 
Critical point are considered it te"fTound’ that /or ort oo 
for hydrogen peroxide. Thus we get lower bounds for 1/6 
in both two and three dimensions and we can use these 
Velues as cOntidence Jimitss TALL the contidences amis 
Quoted in This se€etion are a result of this: mevned or 
Gher ratio mevhod. 

Tables 9.5 and 9.6 are a summary of the results of 
Pusse Varvous Mevnods. he number son Csi eniitesn te letras 
quoted represents the apparent precision of the various 
SonstsvenG Lesulvs and has-novching vo do with tHe factual 
accuracy of the results. -Since many of these seemingly 
very precise results differ from lattice to lattice the 
Various Estimates of the exponents have not converged 
as much as the methods seem to indicate. For this reason 
no confidence limits are quoted where the ratio method 
and method. 4 eave mo results. 

The two dimensional results agree remarkably with 
Sealing onjall paths but the diagonal. ingaierc aces 
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Table 9.4 
Value of (u = 1)(d/du) loz M at the critical moins 
l= Plousing each successive coefficients of Ghe critical 


isotherm magnetization series for the honeycomb and hydro- 


Penyperoxside lattices. 


Degree of Polynomial Numerical value of polynomial 


1 0.05742 0.10951 
2 Preemie 262 0.13546 
3 0.06199 0.15098 
4 0.05884 0.16053 
0.06318 0.16643 

6 OF0b518 0.16997 
‘i 0.06503 ORR Lon 
8 0.06348 On der4 
) 06486 Oda 3oe 
10 .06554 0.17522 
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14 .06561 0.17993 
ae 06561 Oh, Intelelre ye 
Wa .06570 OL uikenisgl 
17 06577 Oy lo2n 7 
18 06584 0.18282 
19 06590 OF es 4a 
20 06590 0.18396 
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Percent of scaling predictions. 

In three dimensions the results also agree quite 
well with scaling but are of an order of magnitude less 
precise. This is due to the imprecision in the estimate 


or the critical point and the critical) exponents 6 and 


Noe 


9.2 Analysis of Diagonal Series 


ihe analysis of the diagonal serves for aly tatircees 
REPresents a very “difitceuly problem. Hieure 94s) 214 
tavio plou Tor the moneycomb latrtace., VAlso'snown gare 
straight lines corresponding to 8, = a late B, = A she 
and Bs = 1/7. In Table 9.7 a sequence of approximates 
LOmLne exponent. civchvasValue 101 CNe erie. cal spony 
i =e as tabulated... (iiom sone Hioure, Sas eand Tap le so. 7 
PMere 1S NO eViIGence Uhat with the number of Vrermnssava lua 
bles che scaling value of, B= V15 holds and the spontaneous 
Maenetization value of 1/73 is also implausible. The ratios 
do however fit rather well to the line corresponding to 
Bo ou lyiy? 

When method 1 is used on this series the results 
are more consistent with scaling theory. In Table 9.8 
the Jocations of the poles and resulting residues are 
tabulated for @ few higher and central Padé approximants 


to the logarithmic derivative of the magnetization. 
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PIGURE 9.3 


RATIOS Ln VS. 1/n FOR THE HONEYCOMB 


MAGNETIZATION ON THE DIAGONAL PATH. 
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Sequence of approximations to the exponent given 


Coe crivircal poinc uae 1 for the honeycomb magnetization 


Chie vie wdleconal path. 


Me) Gy S © Wy 2 & MW f= 


= 
S 


bite 7 


= im 


San Ly 


SeTa.0 


By iL 0 


0070.0 
Sel. 


| ao 
iets entnicee ‘sie B28 oes np Je 
nk ema fone he vOay acorn ia #0" ari 


gy ‘ | 1 
md : ova) Hf j - eee a, 


. 
=— 


en oe 
mgt intot am 


+ ersg L ‘ _ pt 


a "hy a 


83 


Table 9.6 


Padé approximants to Magnetizavion on the diaronal 


pac on the honeycomb lattice. 


Approximant d/du log M : 
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The two estimates are plotted in Figure 9.4 and a very 
Smooth curve results. The intersection of this curve 
Wits = ls fat the point 0.0765. 

A short evaluation table of Padé approximants to 
the function (9.1) at the critical point pw = 1 for the 
honeycomb diagonal series, is given in Table 9.9. 
"best" estimate of the critical exponent from this 
table is 8, = OL0765. 

When method 3 is used on this series the power 
Woveh reproduces: the Known critical point the best is 
=1j.0 Which corresponds vO a Critical exponent or 
0.0769. In Table 9.10 the evaluation of the series 
Co eruneGion (9.1) at the critical pointy =. Using 
Successively higher cCoefiicient of the series is cabu— 
laced or the honeycomb: diaconal series. Neve tne 
Sequence is slowly convergent buv very regular. ihre 
makes the upper bound of 0.07968 very probable. 

When these five methods are used on the derivative 
Govthe diagonal series, the estimave Tor tunes critvical 
exponent of the derivatives is -0.85714 for all five 
methods. This corresponds to a critical exponent for 
the magnetization on the diagonal of 0.14294 + 0.0001 
Or 1/7. The results seem to tall inve two contradictory 
eroups. For the analysis of the derivative and the ra- 
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FIGURE 9.4 


LOCATION OF THE POLE Vs. THE RESIDUE AS 
DETERMINED FROM PADE APPROXIMANTS TO 
(d/du)log M(u) ON THE HONEYCOMB DIAGONAL 
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Table 9.9 


Evaluation of Pade approximants to (l-u)(d/dp) 
loge sMti eat the critical point. u = 1 for the diagonal 


seraes On the honeycomb Jattice, 


Approximant 
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Table 9.10 
Value of (uw- 1)(d/dp) toe MU) at the critical 


point uw = 1 using each successive coefficient of the 


diagonal series for the honeycomb series. 


Deeree sot Polynomial Numerical value of 

evaluated polynomial 
i 0.07695 

2 0.09279 

2 0.09057 

4 On06260 

5 0.08704 

6 O.War29 

7 0.08312 

8 0.08429 

9 0.08376 

16) | 0.08238 

11 0.08233 

We 0.08194 

ie 0.08123 

14 0.08104 

15 0.08062 

aus 0.08026 
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very conclusively. For all analysis on the magnetization 
other than the ratio method, 0.0765 is Ziven Very consise— 
tently and method 4 seems to conclusively rule out By = 
Die ee eo thermodynamics rules out a i. sane] 
17? > 1/8 = 8. This grouping of the results of the 
anaiveis intora 1/7 and a 1/13 is 4 feature cof all the 
other two dimensional lattices studied. The diagonal 
series on the three dimensional lattices have a similar 
DUG tess marked =erouping of the resuliys of che ena iyci.. 
A possible reason for this. contradictory @roupine 
can, be Tound by a “study of all the roots, to the Pade 
approximants to the logarithmic derivative of the honey- 
- comb diagonal series. Padé analysis of the logarithmic 
derivative reveals the pattern of singularitaes alius— 
traved In Figure 955. Note -all the non-physical isineus 
iariti1es are (Of Unev order Cf nic distance Prem eine 
Ori ein, and. this pavvern mignu atlecy tne Taviocrvern, 
seriously. Also the pole just beyond unity om the real 
axis shows a very strange behavior. This pole appears 
to be converging on the point s = 1. Through the five 
highest degree of Padé approximants Cie. OO Lesmevies 
Steadiily from 1.30te 12, while the location of svulc 
physical pole at s = 1 stays fixed. This pole has a 
positive residue that is about an eighth of that of the 


physical pole but is growing rapidly as JLue pole 
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FIGURE 9.5 


SINGULARITIES OF THE DIAGONAL SERIES ON THE HONEYCOMB 
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approaches unity. This behavior seems to indicate the 
presence of two confluent singularities. 

Also when the derivative of the diagonal series 
Pe mult p led egy et 1) = ae to eliminate the apparent 
pole at s\= 1, the analysis shows that there is a Very 
Sunt Cant Simei lariay solu WeMed Mati sera Gece — ae 
This also seems to support the presence of two confluent 
poles. The grouping of the values of the exponent into 
two groups can also be explained by the presence of a 
confluent pole. For example the following function 


might give the behavior 


17 


Mai = eCCn = RAG eye Nee > 


The Patios of Une series expansion to this) tunc— 
Drone wal be dominated, an the Carlier Germs py tne 
exponens =-1/7. Since l// >> 170 = 1715 the contribu= 
tion from Ghe 1/7 singularicvy might vovalily mask the 
IVS singularity an the derivarive buy mot in ehe mdag- 
nevigation itself. The occurrence of two conrluent 
Singularities has been the assumed behavior of vhe 
diagonal series and the value closest to the scaling 
value has been accepted as the "best" value for the 
critical exponent. Evidence has been found, on the 
bathe for p= 35 2.12/72 and72, which indicates the 


occurrence of 2a confluent pole for these paths also, 
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Tt is thought that the appearance of a confluent singu- 
larity might be a general feature of all paths other 


trance Zero Liecld and vcriticalewsovlerm pavis. 
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Suitable conformal transformations might sharpen 
estimates of the critical point and the confidence 
limits for many of the series. A good choice of a trans- 
formation of the series expansion variable will smooth 
out the ratios and make function (9.1) regular. Trans- 


Pormavvons sor. phe Torm 


a all set: 
Ss Sl = “1. 2euser Cageets) 


worked well on the diagonal series for the five lattices. 


As an example the transformation 
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eee pacers sate 


Doe emDech Used OW Pole dtaconal Serres ore tie square 
Jattice. The ratios become flav and give am estimate 

Of one Seventh, for the critical exponent. Also, waen 

Chie pransitormation is Wsed on the Square diaconal series, 
ther eVvaluavlon of mne series: vO function (9.41) at. the 
Critical pointy =) 1 using successively Nigher -coctticients 
of the series becomes regular as shown in Table 10.1. 

From this table it can be seen that an upper bound of 
0.07964 te very probable. his technique might also be 
Heed profitably om Une other Series analysed. 


The paths corresponding to 


z/z. « pe (10.3) 
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Taole ao. 
Neaiue ot (s*-1)(d/ds*)loe MCs*) ac the crivacat 
point u = s*= 1 using each successive coefficient of 
the transformed diagonal series on the square lattice, 


when wp = s = 0.84 s¥*/(1 - 0.16 s#). 


Degree or Polynomial Numerical Value of 
evaluated polynomial 
7 Ces Zaks 
8 0.08957 
g Ogee! 
10 eWay Ges 
teak 0.08348 
L2 0.08481 
is Ome oa 
14 0.08355 
1) 0.08249 
16 0.08252 
Ae yi 0.08190 
18 . OL08l7 6 
19 Onna 2 
20 et Ohe pilelieah 
al, 0.08077 
22 0.08056 
sie) 0,08029 
oy 0.08008 
25 0.07984 
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might also be analysed for small integral values of p 
and 1/p. Near the critical point these curves become 
Straight lines of slope p. Therefore, according to 
Scaling,these paths should have the same exponent as 
the diagonal path. The analysis of the series on these 
paths will give estimates of the critical exponents just 
above and below the diagonal curve. The preliminary 
analysis indicates that these paths ell havera crivical: 
exponent Of 2/15 plus orsminus: 57.) “This 2ives more 
Validity to che Conclusion that. scaling Gheorvyenoltds. on 
Gaewolasonal pata. 

mune Eftecyu of, Coniduenl singularities; On the Tratioc 
and Padé approximants to truncated series has never been 
Scvudved an-any detalii. Thevconclusi ons wot. 1nis* Chess 
show the need for such studies on known functions. 

Ali. the analysis completed so Tar agrees very weld 
With the predictions of scaling theory. Thus one more 
positive test of the validity of the scaling hypothesis 


has been added. 
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APPENDIX 


Ising model low temperature magnetization series 
on paths Yof the form Z/Z = fee) enor phen hconey— 
comb, Square, trianguler,. hydrogen peroxide, and diamond 


lattices. 
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